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ABSTRACT 

Atmospheric mass-loss from Hot-jupiters can be large due to the close proximity of these planets 
to their host star and the strong radiation the planetary atmosphere receives. On Earth, a major 
contribution to the acceleration of atmospheric ions comes from the vertical separation of ions and 
electrons, and the generation of the ambipolar electric field. This process, known as the "polar wind", 
is responsible for the transport of ionospheric constituents to the Earth's magnetosphere, where they 
are well observed. The polar wind can also be enhanced by a relatively small fraction of super-thermal 
electrons (photoelectrons) generated by photoionization. We formulate a simplified calculation of the 
effect of the ambipolar electric field and the photoelectrons on the ion scale-height in a generalized 
manner. We find that the ion scale-height can be increased by a factor of 2-15 due to the polar wind 
effects. We also estimate a lower limit of an order of magnitude increase of the ion density and the 
atmospheric mass-loss rate when polar wind effects are included. 
Subject headings: planets and satellites: atmospheres 



1. INTRODUCTION 

In the past two decades and in particular, following 
the Kepler m ission, hundreds of exoplanets hav e been 
detected (e.g., |Schneider|1995| [Mayor et al.|2003| . Many 
of these planets are gas giants observed at an extremely 
close orbit of less than 0.1 AU from their host star (an or- 
bital period of less than 10 days), and are classified under 
the term "Hot-Jupiters" (HJ). The unexpected close-in 
orbit of HJ has stimulated many science investigations 
regarding their formation, evolu tion, and tidal interac- 
tion (e.g., Papaloizou et al. 2007| and references therein) 



Co- 



ure 



their magnetic interaction with the host star (e.g., 
hen et al. 2010 and references therein), a nd the struc 
and d ynamics of their atmospheres (e.g., [Showman et al 
20081 and references therein). 



In such a close orbit (especially if the star and the 
planet are tidally-locked) , HJ are expected to receive 
extreme ly large amounts of stellar X-ray and EUV ra - 
diation ( |Penz et al.]|2008[ |Cecchi-Pestelhm et al.||2009[ ). 
It has been argued that this high EUV radiation can 
lead to a strong photo-evaporation of the planetary at 



mosphere and high mass loss rates ( Lammer et al.||2003 
Baraffe et al.|2004[|2006| , leading to a less massive plan- 
ets. However, this could not be supported by the ob- 



served mass distribution (Hubbard et al. 2007). Obser 



vations of Lya emission from the HD209458 system have 
suggested that the planet occupies an inflated hydrogen 
corona with outflow velocities of 50 — 100 km s _1 and a 
mass-loss rate of about 10 10 g s _1 . However, there is a 
debate on whether the observations are effected by the 
host star or whether the observed features are of plan- 
etary origin (|Vidal Madjai et al. | |2003[ |Ben-Jaffel|[2007[ 
Vidal-Madjar et al. 20081. A more recent observation of 



the system, as well as of the HD189733 syst em reveled 
a smaller mass-loss rate of about 10 8 g s -1 (Lecavelier 
Pes Etangs et al.||2010[ |Linsky et aL"|2010| ). 
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On the theoretical side, several models for atmospheric 
escape from HJ have been developed in recent years. A 
detailed models for the chemistry, photoionizat ion, and 
aeronomy of HJ were d eveloped by |Yelle[ ( |2004 ) and by 
Garcia Munozl j2007[ ). |Tian et al.| ( |2005| ) ancTjMurray- 



Clay et al.| ( |2009p performed hydrodynamic cal culations 



of the r mally driven atmospheric es cape, and IStone fc| 



Proga| ( |2009j ), |Trammell et al.| ( |2011[ ), and |Adams| ( |2011[ ) 
included the planetary magnetic held geometry, which 
confines the escaping gas to regions of open field lines. 
The models above predict mass-loss rates not higher than 
10 10 g s _1 . Some of the models also included the incom- 
ing stellar wind and found that the planetary outflow is 
ought to be suppressed by the wind. Non of the mod- 
els predicts a sufficiently high mass-loss rate so that the 
planet can be evaporated in a relatively short time-scale. 

In the Earth's upper atmosphere (as well as in other 
planets), there is a well observed physical process which 
plays an im portant role in the ac celeration of ions. The 
polar wind ( Banks fc Holzer|1968 ) is the outflow of plan- 
etary ions along open held lines. The main driver for this 
process is the ambibolar electric field, which is propor- 
tional to the electron pressure gradient. Since electrons 
are more mobile than ions, a charge separation is created 
along the magnetic field direction, leading to an electric 
potential that acts on the ions to retain charge neutral- 
ity. The end result is an acceleration of the ions by this 
electric field so that the ions are dragged by the electrons. 

The electric field applies a force propotional to the neg- 
ative gradient of the electron pressure. Using some sim- 
plifications, the resulting force is approximately equiv- 
alent to half of the gravitational force on the major 
ion species and directed oppositely. Since + is the 
major ion species in Earth's upper ionosphere, the re- 
sult is a supersonic flow of H + and an increase in the 
+ scale height. In addition, photoelectrons, which are 
highly energized electrons due to photoionization (the 
tail of the distribution function), can significantly in- 
crease the electron temperature, leading to an enhance- 
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ment of the ions acceleration (Lemaire 1972). In the 
Earth's upper atmosphere, the velocity of is lower 
than the escape velocity. Nevertheless, + is observed 
to serv e as a significant plasma s ource in the mag neto- 



that: 



sphere |Lennartsson et al.||1981 ). Tarn et al 
Tarn et al. 1(119981) have demonstrated by nun 



H1995| | and 
by numerical simu- 
lation that photoelectrons indeed, can accelerate + and 
H + , while they obtained an unrealistic electro n tempera- 
ture o f 40, 00 K. An additional simulation by Kh azanov| 



et al. (19971 resulted in a more realistic electron tem- 
perature of 16, 000 K. Recent numerical simulations by 
Glocer et al. (2012 ) also included the effects of photoelec- 
trons to look at the global outflow solution and compared 
with in-situ observations. Their simulations showed that 
the polar wind mechanism is responsible to the transport 
of ionospheric H + and + , and that only a small frac- 
tion of photoelectrons can significantly contribute to the 
ion acceleration. 

In this letter, we investigate how the ambipolar electric 
field and the fraction of photoelectrons reduce the grav- 
itational potential, and therefore, increase the ion scale- 
height and the ion density at the top of the atmosphere 
of HJ. We also calculate how the mass-loss rate for H + is 
modified by this effect. Due to the high EUV radiation, 
the fraction of photoelectron in the atmospheres of HJ is 
expected to be higher than in the Earth's case, leading 
to a much greater increase of electron temperature. 

In Section [2j we calculate the change of the effective 
gravity and the ion scale-height due to the ambipolar 
electric field and photoelectrons. We present and dis- 
cuss the results in Section |3j and draw our conclusion in 
Section HJ 

2. MODIFICATION OF THE ION SCALE-HEIGHT 
BY THE AMBIPOLAR ELECTRIC FIELD 

In the derivation bellow, wc follow the standard model 
for the polar wind, but we include the effect of the pho- 
toelectrons on the solution. For a planetary atmosphere 
consiststing of electrons, photoelectrons, and ions, charge 
neutrality requires that: 



n e o + n a0 = n i0 , 



(1) 



where n e Q, n a Q, and Um are the electron, photoelectron, 
and ion number densities at some reference altitude, Tq. 
From Eq. [T] we can define the fraction of the photoelec- 
trons, /3, as n a Q = /3riio, and the fraction of electrons as 
n e0 = (1 - f))mo- 

Our goal here is to calculate how the effective grav- 
ity at the surface is modified when taking into account 
the photoelectrons and the ambipolar electric field, and 
investigate how this modified gravity affects the ion scale- 
height Hi . We will compare Hi with the unchanged scale- 
height Hq which contains the surface gravity g but not 
the ambipolar electric field. 

We begin by assuming a hydrostatic ion density profile: 

n i {z)=r Ho e-(*-^l Hi , 



with the ion scale-height, Hi 



hTi 



(2) 

-, where k is the 

Boltzmann constant, Ti is the ion temperature, rrii is the 
ion mass, and g e ff is the effective gravity. Without the 
effects we study here, g e ff — g. Conservation of the 
photoelectron mass along a magnetic flux-tube requires 



n a0 u a0 Ao = n a u a A, 



(3) 



with u a o and u a being the photoelectrons velocities, and 
Aq , and A being the magnetic flux tube cross-sections at 
the reference altitude and at some altitude, respectively. 
This equation implicitly neglects any scattering of the 
photoelectrons. In a magnetic dipole geometry, the mag- 
netic flux conservation requires that AqBq = AB, with 
B = C/r 3 being the dipole field magnitude as a func- 
tion of radius (C is a constant), and Bo = C/r 3 ) is the 
field magnitude at the reference altitude, rQ. Therefore, 
Aq/A = rjj/r 3 , and we have: 



laO (~) . 



(4) 



assuming u a Q — u a as a lower limit. Using Eq. [2] and HI 
the electron density at altitude z, can now be obtained, 
assuming z = 0, r = R p , and r = R p + z: 



n e (z) = m(z) - n a (z) 



n i0 



-z/Hi 



ml 



(R P + zf 



( 5 ) 

The effective gravity is modified by the ambipolar elec- 
tric field as g e ff = g — -^r, with the ambipolar elec- 



tric fi eld (positive for ions) defined as (Schunk & Nagy 
2004|: 



En 
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Here T e is the electron temperature, and e is the electric 
charge. At the planetary surface, z = and so we obtain: 



E\\(z = 0) 
which yields: 

9eff =9 + 



kT„ 
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T e + (1 - p)Ti 
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(8) 



(9) 



In Eq. |9j g is modified by the ion and electron tempera- 
tures, and by the fraction of photoelectrons. For the case 
of p = and Ti = T e , the well-known reduction of the 



effect ive gravity of the ions by half is obtained (|Gombosi 
20041 



Xs shown by previous models ( [Tam et al. 11 1995 



anov et al.|1997[|Tam et al.|1998|1Clocer et al.|20 



Khaz- 
51), the 



electron temperature is highly iflected by even a very 
small fraction of photoelectrons. In our model here, we 
assume that Tj = 1000A". Despite of the higher ion tem- 
perature expected in HJ, the effect studied here is driven 
by the difference between T e and Tj, so that it should 
scale with the increase in Tj. We scale the electron tem- 
perature with the percentage of photoelectrons and Ti 
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using two different models. One is based on the elec- 
tron temperature distribution at t he top of the Earth's 
atmosphere from Khazanov et al. ( 1997 ) : 



T e (/3) =Ti* 2 6+log ^, 
with 1000 < T e < 16000/T for lO^lO" 4 



(10) 



i) < (3 < 

10~ 2 (1%), and a more modest function with 1000 < T e < 
10000.fr: 

T e ((3) =7W.8 6+log/3 . (11) 

With the above models for T e , the modified gravity and 
scale-height can be calculated as a function of the frac- 
tion of photoelectrons. 

3. RESULTS & DISCUSSION 
3.1. Results 

Figure [Tl shows the effective gravity as a function of 
the photoelectron percentage for the two models for T e . 
For (3 = 0, we obtain g e ff/g = 0.5. In Figure [2l we 
show the electron temperature and the ratio of moaified 
to non-modified scale-height as a function of the fraction 
of photoelectrons, assuming wii = m p , the proton mass. 
Here, we show the solution only for photoelectron per- 
centages of 0.0001 — 1. One can see that if the fraction 
of photoelectrons is even less than 1%, the scale-height 
increases by a factor of 2-15. 

A realistic ion density profile cannot be obtained us- 
ing the simplified calculation we present here. In partic- 
ular, we cannot calculate the density profile of the H + 
ions, since they are expected to attain supersonic speeds. 
Therefore, it is hard to estimate the increase in density at 
the top of the atmosphere and the corresponding increase 
in mass-loss rate. Nevertheless, we can use a hydrostatic 
profile to estimate the ion density change at lower al- 
titudes. In Figure [3j we show the ratio of the hydro- 
static density profiles using the modified and unmodified 
scale- heights, respectively, as a function of the fraction 
of photoelectrons for altitudes of 350 km (~ li?o) and 
1000 km (~ 3H ). The density is increased by a factor 
of 2-3 at 350 km and by a factor of 5-15 at 1000 km. At 
higher latitudes, the hydrostatic solution is probably not 
valid anymore and the ratio in Figure [3] will become too 
large, since the density profile for the unmodified scale- 
height goes to zero faster than the one with the modified 
scale- height. 

3.2. Discussion 

In HJ, the extremely strong radiation is expected to 
increase the fraction of photoelectrons. Therefore, the 
electron temperature should be higher than the ion tem- 
perature, despite of the strong heating at the day side, 
so that the mechanism proposed here should still be sig- 
nificant. The effect should be limited at the night side 
due to the lower ionization rate, and it is not clear yet 
how effective the atmospheric day-night circulation is at 
higher latitudes (where the day-night temperature dif- 
ference is smaller than that at the equator), and at high 
altitudes (where the ion acceleration occurs). 

For a magnetized HJ, the mass- loss is expected to take 
place along the magnetic field lines which are open to the 



stellar wind (as demonstrated by Stone fc Proga 2009 - 
Trammell et al.pOll Adams 2011), and tfiat is exactly 



previously shown that the classical polar wind mecha- 
nism together with the addtion of photoelectrons and 
wave-particle interactions is responsible for the trans- 
port of H + and + out of the Earth's atmosphere. By 
lowering the potential barrier, these processes effectively 
lower the the escape velocity. These processes have also 
been speculated to be important at Jupiter and Saturn 
( |Glocer et aLp007j |Nagy et al.|[l986] ); the major ions m 
the upper atmosphere at these planet are and H + . 
In HJ, the r elative ion abundance s are not known, but 



modeling by 



Garcia Muhoz 



(|2007j) shows that H+ , 



He + , C + , and various ionized hydrocarbons are possibly 
present. The polar wind process should apply to each of 
these planets. Indeed the derivations presented here re- 
flects the basic textbo ok deriviation of the classical polar 
wind ( |Gombosi|2004[ ) to which we have added the effect 
of photoelectons! No other planet specific parameters 
are required. Even neglecting the effect of photoelec- 
trons, the polar wind process by itself could significantly 
increase the ion scale height. 

The relative composition affects the polar wind process 
by changing the parallel electric field. This is because the 
parallel electric field was found to increase with mass. In 
the case of no photoelectrons, if was the major ion 
species (such as at Jupiter or Saturn) than the paral- 
lel electric field would excert an upward force approxi- 
mately equal to one half the gravitation force acting on 
i?3~. In this case the scale height of H$ would increase. 
Lighter consituents such as H + would actually have a 
net upward force resulting in eventual supersonic flow. 
Including photoelectrons, increases the electric field and 
intensifies the effect of the polar wind, possibly resulting 
in a net upward force on heavier species. If H + was the 
major ion, the parallel electric field would be reduced, 
but the effect would still be quite significant. 

The simplified model presented here is insufficient to 
predict the detailed change in the ion density profile, 
but it can predict how the scale height changes. We 
show that this change can reach about a factor of 10 at 
lower altitudes. Therefore, it should also increase the 
mass-loss rate by the same amount assuming the same 
surface area, and without changing the ion velocity at 
the top of the atmosphere. The polar wind is expected 
to further accelerate the ions such that, the ion speed 
should increase as well, so the factor of 10 increase is a 
lower limit. 

In order to perform a more detailed calculation of the 
effect of the polar wind on the mass-loss rate of HJ, a 
more det ailed model is needed, such as the polar wind 
model bylGlocer et al.| (|2007[), |Glocer et al.|(2009D, and 



Glocer et al.| (2012), which is similar to tfiat 



t of |Garcia| 
ambipolar 



where the polar wind process takes place. It has been 



but includes the effect of the ambipolar 
electric field and photoelectrons. The derivation and dis- 
cussion contained in this paper, however, demonstrates 
that the polar wind process plays an important role in 
the mass-loss rate of HJ and should be accounted for in 
models. 

4. CONCLUSIONS 

In this paper, we perform a simplified calculation of 
the effect of the ambipolar electric field and atmospheric 
photoelectrons on the planetary ion scale- height. We 
show that this effect can reduce the effective gravity and 
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therefore, enhance the ion acceleration in the region of 
the planetary atmosphere, where magnetic field lines are 
open. We find that a small fraction of photoelectrons 
(less than 1% of the total electrons) can increase the 
ion scale-height by a factor of 2-15. We calculate the hy- 
drostatic density profiles using the modified scale-heights 
and find that the planetary mass-loss rate should increase 
by an order of magnitude at a minimum, even neglecting 
any increase in the ion velocity due to this the process. 



Since the ion acceleration should be enhanced by the pro- 
cess, we expect the increase in mass-loss rate to be even 
greater. A more comprehensive calculation, however, re- 
quires a more detailed modeling effort. 

We thank an unknown referee for her/his review re- 
port, and Jeremy Drake for his help in preparing this 
manuscript. OC is supported by SI Grand Challenges 
grant number 40510254HH0022. ' 
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Fig. 2. — Electron temperature (top) and the ratio of modified to unmodified ion scale-height (bottom) as a function of the photoelectron 
percentage (0.0001 — 1%) for T e model 1 (solid line) and model 2 (dashed line). 
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Fig. 3. — Ratio of non-modified to modified hydrostatic ion densities for T e model 1 (solid line) and model 2 (dashed line) at z = 350 km 
(top) and at z = 1000 km (bottom) as a function of the photoelectron percentage (0.0001 — 1%). 



